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ABSTRACT 

A generic property of the cuspy simulated virialized halos in cold dark matter cos- 
mogenies is that their concentration is inversely correlated with their mass. This be- 
havior has also been confirmed in observations, although differences in the exact form 
and dispersion of this so-called mass-concentration relationship have been reported. 
Some observational studies of massive halos suggest that they are statistically over- 
concentrated with respect to the expectations of ACDM. Here we investigate the 
impact that various published mass-concentration relationships, both from simula- 
tions and derived from observations, would have on other cosmological observables, in 
particular considering upcoming surveys. We find that an integral measure of lensing 
shear, such as counts of peaks from halos, is very sensitive to the relationship between 
mass and concentration at fixed a%, and the disparity between some reported fits is 
much larger than the impact of uncertainty in cts itself. We also briefly assess the 
impact of baryonic physics on cluster scale observables, using state-of-the-art simula- 
tions, concluding that it is unlikely to give rise to the high concentrations reported for 
some clusters. 

Key Vifords: cosmology - gravitational lensing. 



1 INTRODUCTION 

Galaxy clusters are the most massive bound structures in 
the universe, some with masses in excess of lO^^M©. Many 
observables related to clusters are sensitive cosmological 
probes. The evolution of the abundance of rich clusters has a 
strong dependence on ag,, the a mplitude of mass fluctuations 
on the scale of 8/i~^Mpc (e.g. Fan et al. ( ,1997) ), and the 
covariance of cluster counts is very sensitive to primordial 
non-Gaussianity ( |Oguri| pOOg] ); |Cunha et al.| ( |2010[ )). Fur- 
thermore, our understanding of the interplay between bary- 
onic and dark matter in clusters, and of cluster formation 
and evolution can be improved by comparing observations of 
clusters with those formed in cosmological simulations that 



incorporate baryonic physics (e.g. Mead et al. (2010l). 

Numerical simulations of structure formation in cold 
dark matter (CDM) cosmogonies reveal that virialized dark 
matter halos, with masses spanning many orders of magni- 
tude, possess a cusped universal mass density profile, inde- 
pendent of the initial density fluctuations and cosmological 



parameters (Navarro et al. 19961. The density profile pro- 



gressively steepens from the centre to the virial radius Tvir, 
with an isothermal slope at the scale radius r^. A measure 
of the mass concentration of halos is given by Cvir = '"vir/j's- 
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This concentration parameter shows a good correlation with 
virial mass, both in simulations and in observations, and in 
CDM there is a generic trend that more massive halos are 



less concentrated. Wang & White ( 2009 1 have demonstrated 



that there is also a correlation between mass and concentra- 
tion in virialized halos formed in hot dark matter (HDM) 
cosmogenies, but with more massive halos tending to be 
more concentrated. Given that in CDM and HDM struc- 
ture formation proceeds either via hierarchical aggregation 
or monolithic collapse respectively, their work indicates that 
mergers are not the dominant mechanism responsible for the 
universal profiles and distributions. 

Some analyses of rich galaxy clusters, in particular using 
gravitational lensing data, have shown concentrations higher 



than e xpected in ACDM, e.g. MS 2137-03 (e . g.|Gavazzi et al. 
( [20031 )), CI 0024-^1654 (e.g. [Kneib et al.] (|2()03|)) and~th e 
very well-studied AbeU 1689 (e.g. [Broadhurst et al.|(|2005 1 



Umetsu & Broadhurst (20081; also see Corless et al. (20091 



for a compilation of estimates from various works). Gener- 
ally speaking, factors that may contribute to this apparent 
excess in concentration include having a triaxial halo with 
the major axis close to the line of sight ( |Oguri et al.| ( [20051 ); 



Gavazzi (20051; Corless fc King (20071), where by adopt- 



ing a spherical model (circularly symmetric in projection) 
in the analysis one can bias parameter estimates and un- 
derestimate error bars. Projection of other structures close 



2 



to the line of sight can similarly result in apparently high 
concentration (e.g. King & Corless (2007l). Observational 



results are often compared with dark matter simulations, 
but baryons - with their more complex physics - are dom- 
inant in the centres of galaxy clusters, galaxy groups and 
massive galaxies. One important process due to the presence 
of baryons, is that dark matter halos may undergo contrac- 
tion ([Eggen et al. (19621; Blumenthal et al. (19861; Gnedin 



et al.| ( 2004 1 ) resulting in significant deviation from the form 
arising in dissipationless collapse. 

Differences in the exact form of the so-called mass- 
concentration (hereafter M — c) relationship have been re- 
ported, both measured in CDM simulations and estimated 
from observations (e.g. [Bullock et ah] (|2001[); [Comerford fc| 
Natarajan|(|2007[); [Duffy et al.| (|2008[); |Macci6 et al.| (|2008[); 



Mandelbaum et al.[ ( [2008[ ); [Oguri et al.[ ( [2009[ ); [Okabe et al 
( [2009P ). For example, considering only galaxy clusters whose 
radial profiles are well constrained from combining weak and 
strong lensing data, Oguri et al. (2009) note a 7a excess 
of the concentration of mass compared with expectations 
in ACDM, even including an enhancement for any bias in 



the samples of lensing clusters (e.g. Hcnnawi et al. (20071). 
[Mandelbaum et al.[ ( |2008| performed a statistical analysis to 
estimate the mass-concentration relation over three orders 
of magnitude in mass using the stacked weak lensing sig- 
nals of galaxies, groups and clusters in SDSS. Their analysis 
used weak lensing data in order to avoid the central regions 
of halos where baryons dominate, and to minimise errors 
from misestimation of clusters' centres. They find a M — c 
relationship with a slope consistent with simulations, but an 
amplitude that is about 2cr below theoretical expectations. 

In order to make progress in understanding observed 
departures in the relationship between mass and concen- 
tration - beyond that expected from intrinsic scatter, and 
from known observational biases - in this paper we start by 
considering how particular M — c fits derived from obser- 
vations or measured in CDM simulations would impact on 
other cosmological observables in upcoming surveys. The 
Dark Energy Survey (DE^, for example, is expected to 
contain ~ 3, 000, 000 galaxies out to z ~ 1 and ~ 20, 000 
clusters in excess of 2x lO'^^Mo, thus containing a wealth of 
information on virialized objects. The ingredients of the halo 
model formalism, used in the predictions, are summarised in 
Section 2, in particular the M — c relationships used in the 
later sections. The resulting cosmic shear power spectra and 
the lensing selected halo counts are considered in Sections 
3 and 4 respectively. In Section 5 the impact of baryons on 
cluster-scale mass profiles is explored, in comparison with 
analytic profiles of that mass and with concentrations de- 
rived from M — c relationships. We summarise and discuss 
our findings in Section 6. Throughout we assume a fiat cos- 
mological model, and a dark energy equation of state pa- 
rameter w = —1. 
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THE HALO MODEL DESCRIPTION OF THE 
MATTER POWER SPECTRUM 



In the Press & Schechter model ( Press fc Schechter|[1974 1, 
all matter in the universe is structured in virialized dark 
matter halos. Under the halo model formalism, the halo 
mass function and halo clustering strength, as well as the 
distribution of mass within the halos is used to estimate 



statistical properties of the density field ( McClelland & Silk 



(19771; Scherrer fc Bertschinger (19911; Seljak (20001; Pea- 



cock & Smith ( 2000) ; Scoccimarro et al. (20011; Cooray fc 
Sheth (20021) which readily allows estimation of the im- 



pact of the properties of halos on observables (e.g. Hamana 
eFaLl poMl); INeyrinck fc Szapudi[ (|2007l; [Fedeh & Moscar 



dini (20101; Pielorz et al. (2010l) 



The matter power spectrum P{k), at wavenumber k, is 
the sum of one-halo, P^''(fc), and two-halo, P^'^{k), terms: 



P{k) ^ P^''{k) + P^''{k). 



(1) 



The one-halo term depends on the distribution of mass 
in individual halos and dominates on small scales: 

(2) 

where M denotes halo mass and dn/dM is the halo mass 
function, giving the number of halos in the mass range M — > 
M + dM at a certain redshift, per unit volume, po is the 
present mean mass density of the universe and y{k, M) is 
the Fourier transform of the matter density profile of halos, 
normalised by their mass. 

The two-halo term arises from elements in distinct ha- 
los: 



P^''{k) 



P""(fc) 



dM 



dn 
dM 



\y{k,M)\bn{M) 



(3) 

where P^'°(fc) is the linear matter power spectrum, and 
bh{M) is the halo bias that takes account of massive ha- 
los clustering more strongly than low mass halos. 

2.1 Halo mass function and bias model 

The halo mass function and bias are taken from ISheth fcl 



Tormen 



dn 
dM 



dM 



(19991. The mass function is given by 
Po 



M 
Po 
M' 



f{v)dv 



(4) 



A 



1 + 



1 



where v = 



a{M, z] 



1 2 



exp 



-(ai/)/2 



where Sc{z) is the cosmology-dependent critical value of the 
overdensity required for spherical collapse at z, for which 
we use a fitting function (see appendix of [Henry] ( [2000[ ) and 
references therein) 



Sci^) = ^^^S^ [1 " 0.0123 log(l + X')] 



X 



20 



(5) 



1 + z 

appropriate for Sljf = 0. u{M, z) denotes the variance of the 
linear density field, spherical top-hat smoothed in spheres 
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enclosing mass M, i.e. within radius R{M) = (3M/47rpo)^'^^- 
The coefficients a = 0.707 and p = 0.3 are fit from N-body 
simulations, and A is a normalisation constant obtained by 
the requirement that all mass is in halos 

.AN 



Po 



dMM- 



dM 



(6) 



Halos are biased tracers of the background dark matter field, 
being associated with peaks of the initial density field. Under 



the peak-background split ansatz (Kaiseri (1984); Bardccri 



etaL] ( [T986| ); |Cole fc Kaiser] ( [T989| ); |Mo fc White, ( ,1996, )), 
the large-scale halo bias can be related to the halo mass 
function; for the Sheth-Tormen mass function this reads 



bh{M) 



+ 



2p 



5c [1 + {aiy)p] 



(7) 



2.2 Halo mass density profile 



The NFW profile (Navarro et al. (19961) is taken as the 



fiducial halo matter density profile that enters as an ingre- 
dient in the halo model via its Fourier transform. The mass 
density p as a function of radius r is given by 



p{r) = 



(r/r,)(l + r/rs)2 ' 



(8) 



where ps is the matter density at the scale radius rg. The 
concentration parameter c is given by c = Vvii/vs, where r^ii 
is the virial radius of the halo (roughly equivalent to r^oo, 
the radius within which the mean enclosed density is 200 
times the critical density of the universe at the redshift at 
which the halo virialized). The virial mass of a halo with the 
NFW profile is 

„3 



Mvir = 



log(l 



1 + C 



(9) 



which, from the spherical collapse model, is equivalent to 

Mvir = ^Trr^irSvir{z)po (10) 

with (5vir being the virial overdensity. A fitting function for 



the virial overdensity is (Eke et al. ( 1996 1; Nakamura & Suto 
1997) 



ISvr^ (1.0 0.4093 X^ '^^"^) 



(11) 



where X is defined as in Eq. 5. The mass-normalised Fourier 
transform of a density profile p{r), out to the virial radius 
is given by 



y(fc,M) = - 



drp(r)47r r' 



sin(fcr) 
kr 



(12) 



and for the truncated NFW profile this reads (e.g. |Scocci-| 
marro et ah] ( |2001[ )): 

t/(fc, M) = A-Kporl [cos(fcr,) [Ci(fcr,(l + c)) - Ci(fcr,)] 

+sin(fcre) [Si(fcr4l + c)) - Si(fcrs)] - ^^(iTc) ] ^^^^ 
where Ci and Si are cosine and sine integrals respectively 



Si(2) 



Ci(z) = - 



dt. (14) 



For simplicity and ease of comparison with other work, we 
focus on the NFW profile here, although higher resolution 



simulations of halo formation show a softening of the pro- 
file with a logarithmic slope that decreases with radius, the 



so-called Einasto profile (Stadel et al. 20091. In Section 5 



deviations from the NFW profile due to the inclusion of 
baryons are considered. 

The lensing properties of the non-truncated and trun- 
cated NFW halo are obtained by integrating the density pro- 
file along the line of sight to obtain the surface mass density, 
E, often given in terms of a scaled projected distance from 
the centre x = r' /vs where r' is the projected distance. This 
yields the convergence, k,{x), when scaled by the critical sur- 
face mass density Ecrit that depends on the lens and source 
distances. Since the surface mass density is circularly sym- 
metric, the shear 7 can be obtained using the relationship 
7(3;:) = R{x) — k{x), where R{x) is the mean convergence 
inside x. See for example Bartelmann (1996) and Wright & 
Brainerd (2000) for the non-truncated and Takada & Jain 
(2003) for the truncated case (where the halo is assumed to 
extend to e.g. the virial radius). Of importance here is that 
K and 7 depend on the shape of the profile as well as the 



mass. For example Hamana et al. (20041 note that k scales 
approximately oc c for c < 5, and oc ° for larger c. 

2.3 Mass-concentration relationships 

The relationship between the mass and the concentration 
of virialized halos has been determined by various authors 
using N-body simulations or observational data. Here several 
representative fits from the literature are adopted, working 
in terms of virial concentrations and masses as noted above. 
Where necessary, normalisations are recast relative to a halo 
of lO^''/i~^M0 for ease of comparison. In Fig. [l] the z — 
M^ir — c^ir fits are collated. 

• Using N-body simulations carried out assuming the 
WMAP 5-year cosmological parameters (Komatsu et al. 
^2009) , Duffy et al. (2008) derived the mean concentration of 
their full sample of halos (including relaxed and unrelaxed 
halos) as a function of mass and redshift: 



DUFFY X 
C (M, z) 



5.72 



Mvir 



1014/i-1Mq 



(15) 



(1 + 2)0-71 

where = 0.26, = 0.74, h = 0.72 and as = 0.8. All 
their fits used a pivot mass of Mpivot = 2 x 10^^h~^ Mq 
in order to minimise the covariance between the normali- 
sation and the exponent of the mass dependence. The la 
confidence intervals on normalisation, exponent of mass de- 
pendence and exponent of redshift dependence are roughly 
2%, 7% and 6% of their fit values. 

• Comerford & Natarajan (2007) presented an observed 
M — c relation, drawing values from the literature (both lens- 
ing and X-ray determinations for mass and concentration), 
as well as using newly determined masses and concentra- 
tions for 10 strong lensing clusters. Where necessary they 
converted their compiled values to a flat ACDM cosmology 
with fim = 0.3, Q.A ~ 0.7 with h — 0.7. Using this extensive 
sample of 62 clusters they fit a concentration normalisation, 
and an exponent for the mass dependence, keeping an in- 
verse scaling with (l + z) for their fit and using a pivot mass 
of Mpivot = 1.3 X lO^^ft-^Mo: 

in T r j,r 1 -0.15 

c^°^(M,2) = 



10.7 



A'/vir 



(16) 
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Figure 1. M — c relationships at 2 = as described in the text. 



with an error of ~ 40% on the normahsation and ~ 85% on 
the exponent of the mass dependence. 

• Oguri et al. (2009) considered an observed sample of 
10 lensing clusters whose profiles are well-fit by a combined 
weak and strong lensing analysis. Assuming the WMAP 5- 
year cosmological parameters as above, and taking the mass 
and redshift scaling from the fit to N-body simulations of 
Duffy et al. (2008): 



OGURI 



{M,z) = 



14.55 



Mvii- 
lO"ft-iM0 



(17) 



They note that there is a 7a excess of the concentration pa- 
rameter above the prediction from the ACDM simulations, 
when the ten clusters with weak and strong lensing analysis 
available are combined. 

• Okabe et al. (2010) used Subaru data to carry out weak 
lensing analyses of 30 X-ray selected galaxy clusters from 
the Local Cluster Substructure SurvejPl Throughout their 
analysis they assumed Qm = 0.27, Oa = 0.73, h = 0.72, con- 
sistent with Komatsu et al. (20101. The M — c relationship 



was derived from the 19 spectroscopically confirmed clusters 
with observations in 2 filters, that are well-fit by the NFW 
profile 



.75 



1014/i-iM 



(18) 



The authors note that none of the massive clusters in this 
sample show very high concentrations, as has been reported 
for some clusters. They also point out that the scaling of con- 
centration parameter with mass is tentatively steeper than 
predicted in simulations, though there is an error of « 50% 
on the mass exponent. For the purposes of Fig.[l] we include 
an inverse dependence of c on (1 -I- z) and rescale to z = 0, 
given the mean redshift of their clusters. 



3 WEAK LENSING CONVERGENCE POWER 
SPECTRUM 

The power spectrum of the lensing convergence is obtained 
from a weighted integral of the 3D matter power spectrum 
as described in e.g. Bartelmann & Schneider (2001) 



P4e) = 



Ho 



fKix) 



.X 



w^ix) 



0?{X) ' 
(19) 

where x is the comoving distance, Xh is the comoving dis- 
tance to the horizon, and fxix) is the comoving angular 
diameter distance]^ The scale factor of the universe a{x) is 
normalised to unity today. The wave vector k is related to 
the angular wave vector £ and fK{x) through k = ^//k(x)- 
The 3D matter power spectrum is calculated using the halo 
model formalism as described in Section 2. W{x) is the 
weighting function that takes into account the redshift dis- 
tribution of the sources from which the weak lensing signal 
is being measured, along with the relative geometry of these 
sources and the density fluctuations responsible for lensing: 



W{x) = / dx'Gix) 



IkU-x) 

fKix') 



(20) 



where Q(x') is the normalised source distance distribution. 
Note that the lower integration limit arises since only density 
fluctuations at x < x' contribute to the convergence. 

We consider sources distributed in redshift according to 
a normalised probability distribution (Brainerd et al. 1996) 



p{z) = 



V{2./P)zl 



exp 



(21) 



where zq is related to the depth of the observations and 
P controls the fall-off at high redshift. For example, as in 
Lombardi, Schneider & Morales-Merino (2002) when /3 = 1.5 
and 2o = 0.7, the distribution has z ~ 1.05; this is taken to 
be the fiducial distribution and is shown in Fig.[2] The mean 
redshift is representative of current cosmic shear surveys. 

In Fig.|3] the cosmic shear power spectrum for each of 
the prescriptions for the M — c relationships is plotted, as- 
suming the fiducial redshift distribution. The uncertainty in 
the power spectrum, assuming Gaussian errors, is given by 



Kaiser ( 1998 1 



AP4^) = 



(2^+l)/,ky 



(7L) 



(22) 



where fsky is the fraction of the sky covered by the survey, 
the mean-square intrinsic ellipticity is {7int) ssid the number 
of galaxies per steradian for which an ellipticity can be mea- 
sured is n. The first term is the cosmic variance, dominant 
on large scales, and the second term is Poisson noise that 
dominates on small scales. We take {7int)^''^ = 0.4. 

We now focus on the Duffy et al. (from simulations) and 
Oguri et al. (from observations) power spectra, since they 
assume the same cosmology in their analysis, and are rep- 
resentative of relationships reported in the literature. Fig.|4] 
reproduces the power spectra from Fig.jS] now with error 
bars corresponding to a survey of 5000 sq. degrees and (i) a 



PI: Smith; http://www.sr.bham.ac.uk/Iocuss 



^ Since the spatial curvature of the universe is taken to be ii" = 
it follows that /a'(x) = X- 
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Figure 2. The fiducial redshift distribution corresponding to 
Eq.[2T]with 2 Ri 1.05. 
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Figure 4. The upper and lower panels show the band power di- 
mensionless convergence power spectra for the Duffy et al. (lower 
set of points in both panels) and Oguri et al. (upper set of points 
in both panels) M — c relationships. The upper panel is for a 
source number density of 10 arcmin"'^ and sky coverage of 5000 
sq. degrees. The lower panel is for a source number density of 165 
arcmin"^ and the same sky coverage. The horizontal bars de- 
note the binning (described in the text) used to derive the errors 
denoted by the vertical bars. 



Figure 3. The dimensionless convergence power spectra for the 
Duffy et al. (lower solid line), Comerford & Natarajan (upper 
heavy dashed line), Oguri et al. ( shorter dashed line) and Okabe 
et al. (dotted line) M — c relationships with their respective cos- 
mologies (and using the same erg = 0.8 for the three observational 
relationships). 



number density of 10 arcmin" , characteristic of the upcom- 
ing DES, and (ii) to a futuristic survey with the same area 
but a number density of 165 arcmin"^ representative of deep 
space-based observations (e.g. the Hubble Deep Field as 
used in |Huterer] ( |2002p ) . In practice, band power estimates of 
the power spectrum would be made, for which we have taken 
13 equal logarithmic width bins between Z = 10 and 8 x 10'* 
cis in I Schneider et al.| ( |2002[). Th e error estimate of Eq.[22| 
correspondingly scales as~^l7Al. As noted above, this ex- 
pression assumes Gaussian errors, whereas there is also a 



non-Gaussian contribution arising from mode-coupling due 
to non-linear clustering that becomes important at I ~ 1000 



for a survey of the depth we consider (e.g. Scoccimarro et al. 
(1999l), and from the statistical fluctuations in the num- 
ber of halos of a certain mass which are sampled in a fi- 



nite survey volume (Hu & Kravtsov (2003l). As noted by 



Scoccimarro et al. (19991, the line-of-sight projection rel- 
evant for weak lensing tends to reduce the importance of 
non-Gaussianity compared with what is expected in 3-D. 
Takada & Jain ( 2009 1 estimate a factor of ~ 2 degradation 



in the cumulative signal-to-noise ratio of the power spec- 
trum amplitude. We return to this point in Section 6 when 
the results are discussed. 



4 HALO COUNTS 



Kruse & Schneider ( 1999 1 have shown that the weak lens- 



ing selected number counts of halos is a sensitive probe of 
cosmology, and this was tested by |Reblinsky et aL] ( |1999[ ) us- 
ing ray-tracing through N-body simulations. Other authors 
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Figure 5. The number of peaks detected above a threshold 
S/N = V, for the Duffy et al. and Oguri et al. M — c relationships 
for a source number density of 10 arcmin"^ and sky coverage of 
5000 sq. degrees, assuming sources at 2 = 1 and z = 0.5; the 
key shows the M — c relationship and source redshift, with the 
middle line of each set corresponding to the estimated number of 
detections and Poisson errors indicated by the error bars. 



have further explored various aspects of cluster detection 



using weak lensing (e.g. Hamana et al. (2004); Dietrich & 
|Hartlap| p010| ; [Marian et al.| pOlOt ; see also the overview 
in the Introduction of Marian et al.| ( 20101), and searched 
for mass concentrations in data (e.g. [Wittman et al. ( 2006 1 ; 
[Gavazzi fc Soucail] ( [2007t ; [Schirmer et al.| ( |2007f )~^ 

We now estimate the number counts due to the most 
massive halos, adopting the Duffy et al. and Oguri et al. 
M — c relationships. Again we assume that halos can be 
described by NFW profiles, truncated at their virial radii. 
Selecting halos by integral measures of their gravitational 
lensing shear is sensitive not only to their mass function, but 
to how their mass is distributed, so one expects a difference 
in halo counts that depends on how concentrated a halo is at 
fixed mass. We consider a Gaussian smoothed convergence 
map, that results in halos with peak convergence its after 
smoothing with a filter of scale 9s . We ask how many peaks 
would have a signal-to- noise f above a threshold z/th, where 
f is defined as Ks/cts, and as is the noise inside the same 
filter arising from having a finite number density of galaxies 
(n) from which to measure the lensing signal, and from their 
intrinsic ellipticities (Kaiser & Squires 1993) 



7in 



(23) 



For sources at redshift z, halos with and an assumed 
M—c relationship, smoothing on a scale of ds and a threshold 
of ^-th 



N(v> t-th) 



1 

47r 



Ax 



dM ^U{u{M,c,z) 



(24) 

where dV/dx is the comoving volume element, 'H{x) is the 
Heaviside step function ('H(j/(M, c, z) — i^th) = 1 iff a; > 0) 
through which the dependence on the M — c relation enters. 

We consider a survey of 5000 sq. degrees and a number 
density of 10 galaxies per arcmin^ usable for shear mea- 
surements, with the data Gaussian smoothed on a scale of 
1', appropriate to the scale of massive halos. In practice a 
filter that is better matched to the profile expected could 



Figure 6. The number of peaks detected above a threshold 
S/N = u, for the Duffy et al. and Oguri et al. M — c relationships 
for a source number density of 10 arcmin"^ and sky coverage 
of 5000 sq. degrees, assuming sources at z = 1. For the WMAP 
5-year parameters (trg = 0.8) the lower (upper) shaded regions 
indicate ztlcr offsets in concentration normalisation (from intrin- 
sic scatter) in the Duffy et al. (Oguri et al.) M — c relationships. 
For the Duffy et al. relation, the dotted lines indicate the impact 
of offsets in erg approximately corresponding to the 2cr bounds 
(erg = 0.77,0.83) from Komatsu et al. (2010). 



also be employed as in e.g. Bartelmann et al. (20011, or as 



discussed in Section 6.2 an optimal filter tailored to min- 
imise confusion of real halos with large scale structure, e.g. 



|Maturi et al. ( 2005 1 , but for the purposes of illustrating the 
differences between M — c relationships, Gaussian smoothing 
suffices. We first assume that sources are at each of z = 0.5 
and 2 = 1.0 and fix the cosmology at the WMAP 5-year 
parameters (erg = 0.8) that were used in the Duffy et al. 
simulations, and adopted in the Oguri et al. analysis. The 
results are shown in Fig. [5] 

Assuming a scatter in the M — c relations of CT(logjo c) = 
0.15 as estimated in Duffy et al. (2008), Fig.[6]shows predic- 
tions for ±lcr offsets in the mean relations, now focusing on 
sources at z = 1. Finally, the comparative sensitivity to erg 
is also illustrated in Fig. [6] with predictions for the Duffy et 
al. M — c relationship, and now taking erg = 0.77,0.83, an 
offset roughly consistent with the current 2cr bounds from 
[Komatsu et al.| ( |2010"| . 

We keep sources at 2 = 1, and for ag = 0.8, flm ~ 
0.26, consider the dependence of counts on various of the fit 
parameters, for a general M — c relationship 



c(M, z) = cnorm (1 -I- z) 



M 

1014/i-iMq 



(25) 



Note that this is for the purposes of illustrating how the 
counts would vary for changes in two of the parameters, and 
that some of the combinations are rather extreme. For a 
threshold Vth ~ 4 the counts are shown in Fig. [7] 



5 BARYONIC PHYSICS 

The presence of baryons will cause differences in the distri- 
butions of concentrations measured from dark matter only 
simulations, and derived from observational data. One im- 
portant process is baryonic cooling that can lead to adia- 
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Figure 7. The number of peaks detected above a threshold 
S/N = Uth = 4 for a general M — c relationship described in 
the text, for a source number density of 10 arcmin"'^ and sky 
coverage of 5000 sq. degrees, assuming sources at z = 1. The 
upper, middle and lower panels fix zpower=-0.71, cnorm=5.72 
and mpower=-0.081 respectively, and the labels on the contours 
denote the number of peaks. 



ics constraints from the Sloan Lens ACS Survey (SLACS; 
e.g. Bolton et al. (20061) massive (~L,) early-type galaxies 
- that tend to be central galaxies in a group or cluster, or 
lack more luminous neighbours - on average have isothermal 
density profiles on scales of a few kpc, with a small scatter of 
< 10% (e.g. [Koopmans et al.| (|2009|). This is well inside the 



scale radius of a typical object in their sample, where in the 
absence of baryons the profile is expected to be flatter than 
isothermal. Satellite galaxies - with more luminous compan- 
ion(s) - tend to have marginally steeper mass-density profiles 



(e.g. |Auger| ( |2008[ ); |Treu et al.| ( |2009[ )), which from N-body 
simulations is consistent with tidal truncation by an exter- 



nal potential (Dobke et al. (20071). Schulz et al. (20091 fit 



NFW profiles to stacked weak lensing measurements in the 
outskirts of SDSS early-type galaxies, and find that the mass 
inside the half light radius directly obtained from stellar ve- 
locity dispersion observations far exceeds what is expected 
from the extrapolation of the weak lensing-derived profile, 
even when the contribution of stars is also accounted for. 
This is consistent with modification of dark matter profiles 
by baryons. 

Focusing on cluster scales, since these objects contribute 
to the halo counts discussed in Section 4, we now illustrate 
the differences between having dark matter only and bary- 
onic physics (with cooling and star formation, and with 
and without AGN feedback). We use the total mass pro- 
files of a ~ 3 X 10^* M0 galaxy cluster from the simulations 
of |Puchwein et ah] ( |2008[ ) simulated with different types of 
physics, using the numerical prescription described in |SijackTl 
et al. (20071. The stellar mass fraction of the run with AGN 
feedback is consistent with observations ( [Puchwein et aL\ 
(2010| )). Without AGN feedback the stellar mass fraction 
in the cluster's centre is too high, and the contraction of 
the dark matter is greater than seen in the simulations with 
feedback. Other authors have also recently noted that when 
feedback is included at a level that yields stellar mass frac- 
tions consistent with observations, the central density is less 
enhanced than previously thought. This has been confirmed 



on group scales by Duffy et al. ( 2010 1, where in fact the con- 
centration parameters in the runs with baryons can actually 
be less than in the dark matter runs. 

In Fig.jS] we plot the the mass-normalised FFT, y, of 
the spherically averaged density profile for each of the dark 
matter and baryonic physics runs without and with AGN 
feedback. We also plot the FFTs of analytic NFW profiles 
with the same virial radius, but with concentrations calcu- 
lated using the mean for clusters of that mass from the Duffy 
et al. and Oguri et al. M — c relationships. 



batic contraction of dark matter halos, significantly modify- 



ing their mass distributions (Eggen et al. ( 1962 1; Blumenthal 
^taL] ( p86t ; [Gnedin et aT] |2004| ); |White| ^2004[ )); 

Baryonic physics is likely to be most important for 
lower mass objects (< lO^^M©), due to their shorter cool- 
ing timescale ( [Sllk| ( [T977| ); |Rees fc Ost"riker| (m7\ ; |White| 
& Rees (1978l). A stark illustration of this is provided by 



the failure to reproduce the image separation distribution 
of galaxy-scale gravitational lens systems, when the mass 
function is estimated using Press-Schechter theory (or its ex- 
tensions), normalised to match the local abundance of mas- 
sive clusters (see [Kochanek fc White] ( |2001[ ) and references 
therein). Based on gravitational lensing and stellar dynam- 



6 DISCUSSION AND CONCLUSIONS 

Halos formed in CDM simulations have concentrations that 
are anti-correlated with mass. This M — c relationship has 
been reported for both simulations and observations, with 
some differences in the normalisation and in the evolution 
with redshift. We compiled several M — c relationships from 
the literature, and then focused on two representative fits 
from simulations (Duffy et al. 2008) and observations (Oguri 
et al. 2009) that have the same assumed cosmology. Our aim 
was to examine how, if objects in the universe are described 
by one of these forms, this would manifest in other cosmolog- 
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Figure 8. The mass-normalised Fourier transform of various 
profiles described in the text. In order of intersection with 
k=100/i/Mpc, starting with the uppermost, the curves corre- 
spond to the simulated cluster with baryonic physics excluding 
AGN feedback, the analytic NFW with cO'^URl^ the simulated 
cluster with baryonic physics including AGN feedback, the an- 
aljrtic NFW with c^UFFY ^he simulated cluster with dark 
matter only. 



ical observables measured in future surveys. The predictions 
have been undertaken in the framework of the halo model. 



6.1 Cosmic shear power spectrum 

In Section 3, the cosmic shear power spectra were derived 
for a near-future survey (with approximately the character- 
istics of DES) and a more futuristic space-based survey, for 
halos following the Duffy et al. and Oguri et al. M — c re- 
lationships (see Fig.[4|. Since the assumed cosmology was 
the same for each of these, on large scales the signals are 
identical. Differences in halo properties are manifest above 
I ~ 1000 where, as noted earlier, non-Gaussian errors are 
likely to be important. With the typical ~ 2, or more pes- 
simistic ~ 5 degradation in the cumulative signal-to-noise 



expected (Takada & Jain 20091, the normalisation of the 
M — c relationship would be difficult to directly constrain, 
until deep space-based data become available (lower panel 
of Fig.|4]|. However, as proposed by Zentner et al. (20081, 



self-calibration of cosmic shear tomographic data could be 
used to simultaneously extract information on halo struc- 
ture as well as on cosmological parameters, in the regime 
where I < 3000, until non-Gaussian errors dominate. Given 
that the differences in the normalisations of the two rela- 
tionships considered here are apparent at I ~ 1000, it could 
well be that combining tomography with pushing to higher / 
would allow useful constraints on the M — c relationships for 
even the shallow survey. Zentner et al. note that the concen- 
trations of 10'^^/i~^Mq halos at « ~ 0.2 could be constrained 
to better than 10%. 

6.2 Halo counts 

The numbers of lensing shear-selected mass peaks detected 
above a certain threshold were estimated in Section 4. These 
turn out to be rather disparate for the two M — c relations 



considered, already in the near-future survey as shown in 
Fig. [5] In Fig.|6]we show the results for the peak counts as- 
suming that the normalisations of the mean relations are 
offset by ±lcr, based on the intrinsic dispersion that likely 
reflects formation history. For the Duffy et al. relation, we 
also considered changes in erg from 0.8 to 0.77 and 0.83, at 
a level consistent with the 2(j bound from IKomatsu et al.l 
(2010), so changing the halo mass function. The impact on 
number counts of this departure in erg is much less than the 
disparity between the different M — c relations. Our fiducial 
cosmological parameters are fixed at WMAP 5-year values 
(for consistency with the M — c relations), with errors on erg 
consistent with the current 7- year values of |Komatsu et aL] 
( 2010 1; the errors on parameters will be substantially smaller 
after the Planck mission 

In this work we have employed a Gaussian filter, applied 
to the convergence field. In practice, one would filter the re- 
duced shear, the direct observable from the ellipticities of 
weakly lensed galaxies. Various filters have been developed 



and applied to simulations and to real data. Schneider ( 1996 \ 
and |Schneider et al" (|1998') proposed the aperture mass 
statistic, with a compensated weight function that gives an 
optimized signal-to-noise for approximately isothermal mass 



profiles. Padmanabhan et al. (20031, Schirmer et al. (20041 



and |Schirmer| ( |2004[ ) suggested filter functions that more 



closely follow the NFW profile. Hennawi & Spergel ( 2005 1 



developed a tomographic matched filtering scheme for clus- 
ter detection, that in addition uses redshift information for 
the source galaxies. With and without tomographic informa- 
tion, they consider the performance of the Gaussian filter, of 
the aperture mass statistic - with a commonly used polyno- 
mial filter from the family in Schneider et al. (1998), of the 
Padmanabhan et al. (2003) filter approximately tuned to the 
NFW profile, and of the latter modified by an exponential 
cut-off to minimise the impact of large-scale structure. Us- 
ing cosmological simulations to generate weak lensing data, 
for the case without tomographic information, [Hennawi fc| 



Spergel ( 2005 1 find that at moderate detection significance 
(S/N ~ 3.5 - 4.5) the truncated NFW filter is most effec- 
tive at cluster detection, followed by the Gaussian filter. For 
the highest detection significance, the aperture mass statis- 
tic performs best. For a small range of detection S/N, the 
Gaussian filter performs best. 

Unless the S/N threshold is chosen to be high (~ 4 — 5 
depending on the specifications of the survey), or equiva- 
lently for high mass clusters, the contamination from other 
structure close to the line of sight leads to the detection of 



significant spurious peaks in weak lensing data (e.g. Reblin- 
sky fc Bartelmann|(| 1999');'M etzler et al.|(|200 1); White et al. 



( 2002[ ); |Hamana et al.. ( ,2004j | ; [Hennawi fc Spergel (2005^)). 
Optimal filtering for the detection of clusters, as proposed by 
Maturi et al. (20051, dramatically reduces the contribution 



from spurious peaks arising from the large scale structure 
and from other noise due to the intrinsic elUpticity disper- 
sion and positions of galaxies. More recently, [Maturi et al.| 
(20101 developed an analtyic method suitable for quantify- 



ing the degree to which spurious detections are present in a 
weak lensing survey, for various types of filters. Crucially, by 



see [http://www.rssd.esa.int/SA/PLANCK/docs/Bluebook-l 



ESA-SCI(2005 ll_V2.pdf 
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using N-body simulations, they also confirm that using opti- 
mal filters renders the contribution from spurious peaks due 
to large scale structure low on all filter scales, increasing the 
number of real detections by up to an order of magnitude. 
For the analysis of data from a particular survey, an estimate 
of spurious peaks could be made using the formalism of |Ma-| 
turi et al. ( 2010 1 (Table A.l) , or by the more time consuming 



optimal filtering of synthetic weak lensing data sets derived 
from ray-tracing through cosmological N-body simulations, 
with appropriate observational factors for the survey such as 
the number density of galaxies and their redshift distribu- 
tion from which the shear signal can be estimated. What one 
should also keep in mind is that, as discussed by |Hamana| 
et al. (20041 and Fan et al. (20101, peaks arising due to 



noise are also more likely to appear in the vicinity of real 
clusters. Further, in the context of cosmological parameter 
constraint, which is not the subject of this paper, it has been 



shown by Dietrich fc Hartlap (20101 and Wang et al. ( 2009 1 



that the number and properties of all the peaks detected 
in a weak lensing survey, including spurious peaks due to 
large-scale structure, also contain important information. 

Now to illustrate why differences in concentration would 
impact on the shear-selected peak counts, we consider the 
optimal signal-to-noise of halos as a function of their concen- 
tration, at fixed mass. In the limit of a filter that is exactly 
matched to the halo profile, and neglecting large-scale struc- 
ture noise (so concentrating on the most significant cluster 



mass peaks), the signal-to- noise is given by, e.g. Berge et al, 

(Mol, 



which they show for an NFW profile is 
V oc psrl \/G{c) 



(26) 



(27) 



for a given lens redshift and source redshift (or redshift dis- 
tribution), and number density of source galaxies, with 



G(c) 



0.131 0.375 



-h0.388-5xl0~*c-2.8xl0" 



(28) 



c 

In a given cosmology (fixed critical density, and overdensity 
required for collapse) since 

„3 



Ps oc 



ln(l 



(29) 



and comparing halos with the same mass (and r-2oo), the 
ratio of the signal-to-noise of halos with different concentra- 
tions is 



ln(l + c)- 



ln{l-|-c»)- 



(30) 



which gives the signal-to-noise of a halo, v, compared with a 
fiducial halo, u*. In Fig.|9]this is plotted with respect to the 
concentration of the Oguri et al. lQ^^h~^ Mq halo; note that 
the Duffy et al. halo of the same mass would yield ~ 65% 
of the signal-to-noise due to lower concentration. 

With what accuracy can the parameters of the mass- 
concentration relation potentially be constrained using the 
peak counts, for a DES-like survey? Let us assume that the 
Duffy et al. relation is the true, fiducial, description. Con- 
sidering sources to be at 2 = 1, we estimate the range of 




Figure 9. The signal-to-noise witli wliicli a NFW halo of con- 
centration c is detected, u, relative to a halo of concentration 
c* = 14.55, 1/*, as a function of c. 



parameters that would be consistent with the "observed" 
peak counts, by running different grids of parameters and 
considering the most significant peaks with u > 4.5. To 
approximately account for cosmic variance, the error bars 
in Fig. 5 are doubled, although the statistics of highly sig- 
nificant peaks is dominated by shot noise (Hu & Kravstov 
(2003)). From Fig. 5, it is already clear that the halo counts 
from the fiducial model would be inconsistent with an Oguri 
et al. M — c relation. Given these observational error bars, 
and keeping the mass and redshift dependence of the M — c 
relation fixed, the concentration normalisation can be deter- 
mined to within ~ 7%. Fixing the normalisation and redshift 
dependence, the data would be inconsistent with a flat de- 
pendence on halo mass, but would accommodate a relation 
that is ~ 50% flatter or steeper in Afvir. Finally, for the 
evolution with redshift, fixing the normalisation and mass 
dependence, the observed counts are inconsistent with zero 
evolution of the M — c relation with redshift, but consistent 
with an evolution that is « 30% flatter or steeper in (l + z). 

In practice one would explore which method of peak 
detection was best suited to the data set, either by appli- 
cation to simulations, or by testing the recovery of known 
lower mass clusters. Peak counts looks promising as a tech- 
nique that complements the determination of the M — c re- 
lation from detailed studies of individual clusters (e.g. Com- 
erford & Natarajan (2007); Oguri et al. (2009); Okabe et al. 
(2010)), or from statistical studies (e.g. Mandelbaum et al. 
2008), in particular when the uncertainty on cosmological 
parameters would be further reduced by Planck and sup- 
porting data. 



6.3 The recent extension to the halo model by 
Giocoh et al. 



Recently, Giocoli et al. (20101 extended the halo model 



formalism to include realistic substructure within individ- 
ual haloes, and to include scatter in halo concentration at 
fixed mass. With their prescription for substructure, as mo- 



tivated by the findings of N-body simulations (e.g. Gao 
|et al.| ( [2004[ )), the predicted power on very small scales is in- 
creased, bringing the matter power spectrum derived from 
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the halo model formalism into better agreement with cos- 
mological N-body simulations. 

The increase in the contrast between true halo peaks 
and large-scale structure from this extension would likely 
lead to an increase in halo counts compared with standard 
implementations of the halo model. The extension to the 
halo model formalism of Giocoli et al. (2010), may be espe- 
cially pertinent to properly accounting for the weak lensing 
signal due to low-mass haloes, as they note, and predict- 
ing the counts of low-mass haloes. For our comparison of 
the difference between peak counts for the Duffy et al. and 
Oguri et al. M — c relations, any linear structure would be 
the same in both cases, and the higher normalisation M — c 
relation yields more peak detections. This could then be 
further enhanced by the inclusion of substructure. For the 
detection of the highest significance peaks, where the weak 
lensing signal is dominated by massive cluster halos, the con- 
trast against the large-scale structure is less of an issue. In 
addition, simulations show that substructure itself changes 
the weak lensing signal of these massive haloes by only a few 
percent (King et al. (2001 1). Assessing the impact of this ex- 



Indeed, as noted earlier, it has recently been shown that 
the covariance of cluster counts is very sensitive to primor- 



tended halo model formalism on weak lensing peak counts 
is beyond the scope of this paper, and will be the subject 
of future work. At that point, estimating the accuracy with 
which the parameters of the mass-concentration relationship 
can be estimated could be undertaken with a Fisher matrix 
analysis, also incorporating the (photometric) redshift dis- 
tribution of sources from which the weak lensing signal is 
being measured, and combining information from the mea- 
surement of the cosmic shear power-spectrum for example. 

As discussed further by !&^ocoli (2010) their extended 
halo model with the inclusion of substructure is of much im- 
portance in predicting the flexion signal of haloes, and in the 
investigation of strong-lensing properties such as universal 



magnification invariants. Giocoli (20101 also note that their 



formalism could be further refined, for example to include 
the impact of nonlinear bias and to include halo shapes. 

6.4 Non-Gaussianity 

Some degree of non-Gaussianity in the primordial density 
field can arise in even the simplest inflationary models as dis- 
cussed in e.g. |Liddle fc Lytli| ( |2000| , and we now briefly com- 
pare this with the impact of ag. Commonly non-Gaussianity 
is quantified in terms of the parameter /nl such that the 
Bardeen potential $ can be expressed as the sum of a 
Gaussian random field <E>l and a quadratic departure term: 
<E> = <I>L -I- /nlI^l — The tightest constraints cur- 

rently come from CMB data - using WMAP 7-year data. 



Komatsu et al. (20101 find 13 < /nl < 97. Grossi et al 



( 2007 1 carried out dark matter simulations with various val- 
ues of /nl, and found that only modest departures from the 
halo mass function arise even for large positive or negative 
values of /nl. In the context of weak lensing statistics, and 
using the same set of simulations, [Pace et al.| ( |2010[ ) show 
that non-Gaussianity also has a minimal impact: for values 
of /nl within the current bounds, the impact on the cosmic 
shear power spectrum is well within the uncertainty due to 



the current bounds on ug. Pace et al. (20101 also remark 



that since non-Gaussianity is comparatively much more ef- 
fective at small scales, the scale-dependence of observables 
can be used to aid parameter constraint in future surveys. 



dial non-Gaussianity ( |Oguri| ( |2009| ; |Cunha et al] ( |2010| )). 
Although beyond the scope of this paper, it would be inter- 
esting to examine the extent to which the uncertainty on the 
M — c relation could be reduced by self-calibration of shear- 
selected peak counts (similar to the scheme of Zentner et al. 
(2008)) binned in redshift, and allowing for non-Gaussianity 
at a level consistent with forthcoming Planck data. 



6.5 Baryonic physics 

The impact of baryons on dark matter halos as seen in sim- 
ulations, compared with the change in profiles due to chang- 
ing the halo concentration to be consistent with M — c rela- 
tionships at the same mass, was considered in Section 5. In 
Mead et al. (2010) we demonstrated that simulated clusters 
with AGN feedback are less efficient strong lenses (specifi- 
cally in their production of giant arcs) than clusters simu- 
lated without AGN feedback, and are in fact similar to dark 
matter only halos. Here we considered the ~ 3 x 10^* Mq 



cluster from the simulations of Puchwein et al. (20081, also 



used in Mead et al. (2010) where the stellar mass fraction 
is in very good agreement with observations, in the case 
where AGN feedback is accounted for alongside other bary- 
onic processes. As shown in Fig.jS] the corresponding FFTs 
of these density profiles that enter into the halo model pre- 
dictions illustrate that the simulated model with AGN feed- 
back has much less power on small scales than either the 
simulated model without AGN feedback (i.e. where there is 
over-cooling), or of the analytic model with a concentration 
at the mean of the Oguri et al. M — c relationship. From 
this illustration, and from the general findings of Mead et 
al. (2010) and Puchwein et al. (2010) we emphasize that for 
the observational measurements of the M — c relation made 
on cluster scales, or correspondingly for individual clusters 
with very high concentrations, baryons are unlikely to be a 
significant culprit in yielding results that greatly differ from 
the expectations of ACDM. Note that the number of simu- 
lated clusters, groups and galaxies currently available is too 
small to specify a global y{k, M) form, but an interim pos- 
sibility that we defer to later work would be to calibrate 



an analytic form, such as that used in Guillet et al. (20101, 



using a fairly large set of high resolution simulations that re- 
produce the density profiles and stellar mass fractions from 
observations. 



6.6 Conclusions 

Understanding the physical processes, or the observational 
selection effects, that give rise to some clusters having seem- 
ingly high concentrations, and to tension m M — c relations 
derived from various types of observations, will be possible in 
the near future. In particular, we have shown that the counts 
of halos obtained using weak lensing observations, and the 
cosmic shear power spectrum are sensitive to the M — c re- 
lation. Alongside the exciting preparations for observational 
surveys such as DES, much progress is being made in devel- 
oping tools that are very useful in better understanding the 
implications of observations. For example, during the work 
for this paper, Giocoli et al. (2010) presented an extension 
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to the halo model formalism that incorporates substructure 
in halos, as well as allowing for scatter in halo concentration 
at fixed mass. 

It will be very interesting to compare constraints on 
mass profiles and their concentrations from measurements of 
both weak and strong lensing by individual selected clusters 
from upcoming surveys, from carrying out statistical studies 
in DES such as that in SDSS of Mandelbaum et al. (2008) , to 
explore using tomographic cosmic shear with self-calibration 
(Zentner ot al. 2008) and from measurements of weak lensing 
peak counts. 
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